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Abstract 

Stripe phases are observed experimentally in several copper-based high-Tc super- 
conductors near 1/8 hole doping. However, the specific characteristics may vary 
depending on the degree of dopant disorder and the presence or absence of a low- 
temperature tetragonal phase. On the basis of a Hartree-Fock decoupling scheme for 
the t-J model we discuss the diverse behavior of stripe phases. In particular the effect 
of inhomogeneities is investigated in two distinctly different parameter regimes which 
are characterized by the strength of the interaction. We observe that small concen- 
trations of impurities or vortices pin the unidirectional density waves, and dopant 
disorder is capable to stabilize a stripe phase in parameter regimes where homoge- 
neous phases are typically favored in clean systems. The momentum-space results 
exhibit universal features for all coexisting density-wave solutions, nearly unchanged 
even in strongly disordered systems. These coexisting solutions feature generically a 
full energy gap and a particle-hole asymmetry in the density of states. 



I. INTRODUCTION 



Stripe ordering phenomena on the nanoscale seem to be an inherent consequence of 
electronic correlations in high-Tc superconducting materials. They are most prominent close 
to X = 1/8 hole doping [14] but they were also identified within a broader doping range in 
the pseudo-gap regime [5-7]. However, the nature of the stripe order varies significantly for 
different cuprate materials. 

Unidirectional charge- (CDW) and spin-density waves (SDW) have been detected in many 
cuprates by neutron scattering and x-ray experiments [1, 3, 4, 8-17]. However, the details 
are strongly material dependent. Neutron scattering experiments on La2-a;-yNd^Sra;Cu04 
(LNSCO) at X = 1/8 [1, 12-14] revealed static antiferromagnetic (AF) spin-density wave 
order with a period of eight lattice constants and a concomitant charge-density wave with 
half this period. A similar spin structure was found in La2-a^Ba^Cu04 (LBCO) [3, 9, 10, 
15], in La2-a;Ba^Sra;_^Cu04 (LBSCO) (with y — 0.075) [4], and in La2-x-^Eu^Sra;Cu04 
[16, 17], where SDW and CDW coexist at and near x = 1/8. A common feature of these 
cuprates is an anisotropic low-temperature tetragonal (LTT) phase and, in addition, a strong 
dopant disorder; both, the LTT structure and the dopant disorder are supposed to pin 
stripes. However, not all cuprates exhibit a LTT phase or dopant disorder. While every 
chemical doping of the cuprate parent compounds introduces disorder, its impact on the 
superconducting properties depend decisively on the distance between the Cu02 planes 
and the dopants. In La2-a:Sra;Cu04 (LSCO) dopants are randomly positioned close to the 
Cu02 planes generating effective disordered impurity potentials to the in-plane electrons. In 
contrast, oxygen dopants order in CuO chains in YBa2Cu307_(5 (YBCO) that are separated 
from the Cu02 planes by a BaO plane and they are roughly at twice the distance from 
the Cu02 planes than Sr is in LSCO. Hence YBCO is minimally affected by the dopants' 
impurity potentials and is therefore considered as the cleanest material in the cuprate family. 
In pure LSCO, which exhibits no LTT phase, spin stripes have been detected below but 
no CDW order [18, 19]. Although no static stripes have so far been reported for YBCO, 
electron-nematic order is inferred from anisotropics [20] found in neutron scattering [21] and 
thermoelectric transport [22] measurements. In addition, incommensurate spin fluctuations 
are detected in superconducting YBCO and a static CDW appears in the presence of an 
external magnetic fleld [23, 24]. The different behavior of YBCO may arise from both the 
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absence of dopant disorder and the LTT structure. Stripe phenomena in cuprates must 
therefore be considered strongly material specific. Here we present a comprehensive analysis 
of their delicate response to the presence of disorder. 

In the past few years the theoretical understanding of stripe formation in cuprates has 
advanced. Early on neutron-scattering data on LNSCO at hole doping x = 1/8 [1] sug- 
gested the formation of a spin-ladder structure in the Cu02 planes. This structure is built 
from half-filled, three-legged spin ladders, separated by quarter filled chains. The AF spin 
structure changes sign from one ladder to the next, resulting in a wavelength of eight lattice 
constants. This same spin structure was later detected also in LBCO [8]. Indeed, theo- 
retical models based on coupled spin-ladders describe [25, 26] with some success inelastic 
neutron scattering data on LBCO [3]. Already before the experimental discovery of stripes, 
Zaanen and Gunnarsson [24] and Machida [27] predicted the formation of spin stripes in 
doped antiferromagnets from mean-field analyses of the Hubbard model. The spin stripes 
suggested by Tranquada and coworkers were found also in various numerical calculations 
using the Hubbard model [28], the t-J model [29-36] or the spin-fermion model [37], and 
their existence is by now well established. 

A more delicate problem is the understanding of the coexistence of or competition be- 
tween spin- and charge-stripe order and superconductivity. In disordered systems, d-wave 
superconductivity and antiferromagnetism can coexist [38]. On the other hand, within a 
mean-field treatment of the t-J model and variants of it, the above described spin-ladder 
state was found to coexist with a striped form of rf-wave superconductivity [39-42]. This 
superconducting (SC) state is modulated in space with the same period as the spin struc- 
ture; its SC order parameter is minimal in the center of the AF spin ladder and maximal 
in between the spin stripes. Striped superconductivity in this context corresponds to a 
unidirectional pair-density wave (PDW) state [43, 44]. Notably this PDW oscillates with 
twice the wavelength of the accompanying CDW, which is caused by periodic zero-crossings. 
Below we discuss and contrast the PDW with the periodically modulated rf-wave supercon- 
ductivity (mdSC), which lacks a sign change, thus oscillating with the same wavelength as 
the CDW. 

While the emerging spin order in these models is typically the same in a wide range 
of parameters for constant hole density, the stability of coexisting spin order and super- 
conductivity varies strongly. The t-J model displays two distinct limits which have been 
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assessed within a mean- field approach in Ref. [42] : (z) the strong J hmit, referred to as 
the so-called "y-model", where V parametrizes the attractive nearest-neighbor interaction 
and represents the dominating contribution to the interaction (see Sec. II B) and {ii) the 
weak J limit. The latter is well described by a BCS model which includes a repulsive on-site 
interaction and is appropriately named the "[/-model" (see Sec. II A). In both limits, spin 
stripes form. While in the IZ-model the spin ladders are separated by quasi-one dimensional 
(ID) [41], mutually uncorrelated SC stripes, antiferromagnetism is weak in the [/-model and 
coexists with two-dimensional (2D) superconductivity [42]. The cuprates are estimated to 
be placed in between the two limits. Some cuprate materials show characteristics which are 
qualitatively explained by the y-model, whereas others come closer to the physics of the U- 
model. As pointed out in Ref. [42], for clean systems there is a balance between AF and SC 
correlations in a certain hole-doping range within which specific details matter. This regime 
reacts sensitively to dopant disorder and is likely to describe the physics of 214-compounds. 

In this paper we discuss the impact of inhomogeneities on striped superconductors and 
compare the theoretical results with the experimental observations. In Sec. II we introduce 
the U- and the ^/-models as derived in Ref. [42] and recollect their respective mean- field 
theories. In Sec. Ill we discuss the solutions of the [/-model in the presence of disorder 
and vortices, and the effects of disorder on the solutions of the V-model. In Sec. IV we 
summarize and discuss our results. 

II. HAMILTONIAN 

We follow the general idea that the one-band Hubbard model describes well the low-energy 
physics of the Cu02 planes of the cuprates, including antiferromagnetism and superconduc- 
tivity [45]. At strong coupling, a unitary transformation maps the Hubbard model onto 
the t-J model with an AF exchange coupling J — At^ /U [46]. The non-local interaction in 
the t-J model accounts for both superconductivity and antiferromagnetism already on the 
mean-field level [47-51]. Here we use an ansatz introduced by Kagan and Rice, in which 
the projection to exclude doubly occupied sites is replaced by an on-site repulsion term 
[//2^. leading back to the original t-J model in the limit U ^ oo [50]. On 
mean-field level, this ansatz is equivalent to the fully decoupled BCS Hamiltonian Huv with 
attractive nearest-neighbor interaction of strength V and an on-site repulsion of strength 
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[/, if V is identified with J [42]. We employ this model on a square lattice including ran- 
domly positioned on-site impurity potentials V^™^ and a perpendicular orbital magnetic 
field. Specifically the model Hamiltonian reads 



matrix elements between nearest and next-nearest neighbor sites are denoted by Uj — t 
and tij — t' — — 0.4t, respectively. An electron moving in the external magnetic field 
from site j to i acquires the Peierls phase Lpij = (tt/^q) /J*' A(r) • dr, where $o = hc/2e and 
A(r) = (0, xB) is the vector potential in the Landau gauge. The attractive nearest-neighbor 
interaction is parametrized by 1/ > and the chemical potential /i is adjusted to fix the 
electron density n — /N — 1 — where x is the hole concentration 

A mean-field decoupling of the interaction term —V/2 ^^Ziij) a ^la^l-a^j-a^ia^ leads to the 
standard term for BCS type superconductivity, plus an AF interaction term of the form 
— y ^^^^.^ ^(n^^cr)^i,-cr- The latter results in an energy gain V for two electrons on nearest- 
neighbor sites with antiparallel aligned spins. Consequently, the system exhibits an AF 
phase above a critical V (V t) [42]. In this regime double occupancies become very rare. 
Hence the term U/2'^.^ni(jni-a- is small compared to the TZ-term and has no qualitative 
consequences in this scenario. This case is represented by the TZ-model. On the other hand, 
if y <C no AF order originates from V and antiferromagnetism is controlled by U. This 
case is adequately depicted by the U -model [42] . 

A. [/-model 

As discussed above, the interaction —V/2 X^^^^) ^ cj^cj_^c^_^q^ is not important for mag- 
netism in the limit T/ <C but leads to the standard BCS expression for the SC order 
parameter. Thus the Hartree-Fock decouphng of Huv can be restricted to the t/-model: 




T-Lu — 





(2) 
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where erf = {rii^ — ni])/2. In the following we will focus mainly on hole densities at and near 
X = 1/8. The non-local pairing amplitude is defined as 

^ij = -^(cj^c^t)- (3) 

The [/-model has been the starting point for numerous investigations on disorder- [38, 52- 
56] and field-induced antiferromagnetism [38, 56-58]. In the following we fix the pairing 
interaction V — l.bt and set the strength of the non-magnetic impurity potential V^^"^^ 
to 0.9 1. All fields, i.e. A^j, the local electron density (n^), and the local magnetization (erf) 
are calculated self-consistently from the solutions of the associated Bogoliubov-de Gennes 
(BdG) equations. A detailed derivation of the BdG equations is presented in Refs. [38, 58]. 
The rf-wave order parameter on a lattice site i is defined as 

= - (A^.^^ + A^. .-A^.^^-A^. 0, (4) 

I ^ \ i^i-\-x ' i^i—x i,i-\-y 1,1— y J ' V / 

where Afj — A^^e"^"^^^ . In order to describe a possibly appearing modulation of the su- 
perconducting order parameter, as in the pair-density wave state [41, 42, 44, 59, 60], we 
subdivide the pairing order parameter into three contributions [44, 60] 

Af = Ao + Aq e'^-'^ + A_q e"^^'^^ , (5) 

where the wavevector q parametrizes the modulation of the pairing amplitude and reflects 
a finite center-of-mass momentum hq^ of the electron pairs. Aq accounts for a homogeneous 
component of the order parameter. Two distinct types of SC states are possible solutions for 
T-Lu' If Aq = 0, the SC order parameter shows a stripe pattern modulated with wave vector 
q, where the superconducting stripes are separated by channels of zero pairing amplitude 
across which the SC order parameter changes sign. If the absolute values of |Aq| and |A_q| 
are smaller than the homogeneous component |Ao|, Af is always finite and modulated with 
wave vector 2q. Although the "pure" PDW with Aq = was phenomenologically suggested 
to be the ground state of striped high-Tc superconductors [61], all calculations of mean-field 
type so far led to ground states with dominating uniform component Aq, c.f. Refs. [40, 42]. 
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B. V-model 

For large attractive interactions V ^ t the Hartree-Fock decoupled Hamiltonian of Eq. (1) 
is effectively described by the IZ-model Hamiltonian 

-VJ2 McUc,-. + E K"^^cLc,,, (6) 

where A^^ is defined as in Eq. (3). The term (n^cr)cj_^Cj_o- is responsible for antiferromag- 
netism; it is controlled by the same interaction parameter V as the SC order parameter. 
Since in this limit AF correlations are strong enough to separate the system locally into 
AF regions close to half filling and into hole-rich regions far from half filling, the [/-term 
contributes little and is therefore omitted from l-Lv- 

For hole doping around x = 1/8, the solutions of the BdG equations for Hy in the 
impurity- free case are very similar to those of the [/-model [42]. The main difference is 
the stronger magnetization of the magnetic stripes in the IZ-model. Within the AF stripes, 
superconductivity vanishes, and the regions with dominant superconductivity are spatially 
separated into uncorrelated, quasi ID filaments. Therefore the phase relation between them 
is not fixed and the "pure" PDW state is degenerate to a solution with finite Aq. Because 
of the one-dimensional character of the solutions of the l/-model, they barely adjust to 
impurities in comparison to the solutions of the [/-model. This characteristic difference is 
discussed in Sec. IIIB. 



C. Momentum-space quantities 

For a thorough analysis of our results we also discuss the following momentum space 
quantities. With the Fourier transform of the factorized square of the magnetization 

we approximate the magnetic structure factor. In addition we discuss the absolute squares 
of the Fourier transformed d-wave order parameter 
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and the deviation from the average electron density n 

|l5]((n,)-n)e-^^-^'. (9) 



|n(q)p 



In the strong disorder hmit, the suppression of the electron density at the random impurity 
sites dominates the charge distribution even for small impurity potentials. Coherent charge 
oscillations with small amplitudes are therefore barely visible in the Fourier transform n(q). 
Since we are especially interested in these coherent charge modulations we exclude the 
impurity sites in the presence of strong disorder from the sum over sites in the Fourier 
transformation. 

The momentum distribution is given by 

n(k) = J](4ckJ, (10) 

where the fermionic operators c^^ are the Fourier transform of the real-space operators 
q^. The spectral density at the Fermi energy A(k) = A(k, a; = 0) is determined from the 
imaginary part of the retarded Green's function 

A(k,a;) = -^ 5]lm((c,,;ck.)r. (11) 

a 

A(k) allows for the determination of partially reconstructed Fermi surfaces which are charac- 
teristic for inhomogeneous superconductors. Homogeneous rf-wave superconductors exhibit 
a finite gap away from the nodes on the Brillouin zone diagonals. n(k) is therefore continuous 
everywhere except across the nodes. However, stripe formation leads to a reconstruction of 
the Fermi surface attributed to finite pairing-momentum. In this case the pair density P(k) 
(Eq. (12)) is readjusted and, as a result, the superconducting energy gap no longer covers all 
sectors of the underlying Fermi surface. Impurities, on the other hand, induce bound states 
within the superconducting gap, leading to a redistribution of spectral weight into the gap 
which can be observed in A(k). In this process the transition from partially filled to empty 
states in n(k) narrows, implicating a partial reconstruction of the Fermi surface. 
The momentum distribution of the pair density P(k) is defined as [42, 62] 

p^(k)= j]i(c_k+ci,ckt)r (12) 

q 

It measures the average correlation of the two occupied electron states |k, t) and | — k — q, ^) 
and thereby the parts of the Brillouin zone where electron pairing occurs. Eventually, we 



discuss the information contained in 

P5(k) = ^^(a4^q,Ck,) (13) 

qo- 

on the spin state. This quantity is maximal in those parts of the Brillouin zone where 
antiferromagnetism dominates. Thus the competition between superconductivity and anti- 
ferromagnetism can be analyzed by comparing P(k) and ^^(k). 

III. RESULTS 
A. [/-model 

The Hamiltonian (2) gives rise to unidirectional stripes which have been identified in the 
real-space quantities of clean rf-wave superconductors above a critical on-site repulsion Uc 
[42, 52, 58, 63] . For a clean (T/™p = 0) rf-wave superconductor (dSC) the free energy of 
stripe solutions is typically close to the homogeneous dSC solutions. Depending sensitively 
on the initial conditions of the self-consistent BdG calculations we find either stripes or 
homogeneous AF order, both coexisting with superconductivity. In clean systems, regular 
stripe solutions emerge only if the self-consistency loop is started from a striped initial state. 
However, in the presence of perturbations, such as impurities or vortices, or by a rectangular 
lattice geometry, stripes are pinned and the stripe solutions become robust against changes 
of the initial conditions. This relates to the notion that in the unperturbed superconductor 
fiuctuating stripes are present which become static by the pinning to defects [52, 58, 63]. 

In the t/-model below a critical pairing interaction an emerging stripe state exhibits a 
dominant uniform component of the superconducting order parameter, i. e. Aq 7^ in 
Eq. (5), in addition to the finite-momentum pairing amplitudes. In this case the supercon- 
ducting order parameter Af does not feature a sign change and we call this state a modulated 
rf-wave supercondcutor (mdSC). In Ref. [64] it was shown within a momentum space for- 
mulation that a translation-invariant Hamiltonian in zero magnetic field can nevertheless 
support a "pure" PDW groundstate with Aq = for a similar set of parameters but with a 
stronger pairing interaction V > 2.2 1. There, an analytic approximation to Gor'kov's equa- 
tions was considered. Solutions of this kind are also found in the impurity-free real-space 
model used here, although the PDW solution converges into a local energy minimum and is 
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Figure 1. Density modulations of a d-wave superconductor. SC order parameter for 
an (a) impurity-free and an (b) impurity-pinned (1/™? = lOt) stripe solution, (c) Vertical cuts of 
the electron density and the staggered magnetization in the impurity- free case = 0). These 

calculations were performed on a 24a x 24a lattice where a is the lattice constant. Parameters were 
fixed to T = 0.025 1, t' = -0.4 1/ = 1.5 1, x = 1/8, (7 = 3.3t. 

slightly higher in energy then the mdSC solution. Stripe solutions are indeed observed in 
a broad hole-doping range. In the strong disorder limit, however, we find stripe states only 
close to X = 1/8. 

1. Impurity-free stripe solutions 

Choosing sinus-shaped initial values with a wavelength of 8a for the self-consistent fields 
A^j and = (erf), the impurity-free system (V™^ = 0) exhibits horizontal (or vertical) 
stripes in the real-space quantities over a wide hole-doping range. In Fig. 1 the order 
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parameter, the electron density, and the staggered magnetization are shown for x = 1/8. 
Obviously, a CDW, a SDW, and a modulated pair-density (mdSC) (Aq ^ 0) coexist in 
such a striped superconductor. The superconducting order parameter is jinite everywhere, 
modulated only with an amplitude of a few percent of the average order parameter A^. The 
CDW oscillates around the average density n with an amplitude of only 1%. The staggered 
magnetization exhibits the strongest modulation with an amplitude of about 0.5 /iB including 
a periodic sign change in the vertical direction corresponding to anti-phase domain walls 
between the antiferromagnetically ordered stripes. This solution is the groundstate of the 
f/-model for the chosen parameters and initial conditions. 

The periodic modulation of the magnetization in the vertical direction translates into two 
distinct peaks in the magnetic structure factor at the wavevectors = (27r/a) (l/2,l/2±e) 
with e = 1/8. These magnetic ordering wavevectors correspond to a SDW with period 
8a perpendicular to the stripes as is also obvious from the real space pattern of the staggered 
magnetization in Fig. 1 (c). 

The concomitantly emerging CDW (Fig. 1 (c)) and mdSC (Fig. 1 (a)) modulate with half 
the wavelength compared to the fluctuating SDW. Their Fourier transforms peak at qc/p = 
27r/a(0,±(^) with b = 1/4 indicating a vertical oscillation with wavelength AcDW/mdSC — 
4a. The period doubling of the SDW results from the 7r-phase shift that the spin order 
experiences between neighboring charge stripes. In contrast to the PDW which contains 
no uniform q = pairing component, the SC order parameter is here always of equal sign. 
These wavelengths stay fixed in a hole doping range from x = 0.1uptox = 0.15. For 
X > 1/8, the additional holes collect in the already present stripes which deepens the hole 
rich channels thus increasing the oscillation amplitude of the CDW. 

For a further characterization of the impurity-free case we present in Fig. 2 momentum 
space quantities. The momentum distribution n(k) is shown in Fig. 2 (a) and exhibits an 
anisotropy in k^^ and ky. In the presence of horizontal spin stripes, states with wavevectors 
near k ^ (±7r, 0) are less probably occupied than states near k ^ (0, ±7r). 

In a homogeneous c/-wave superconductor, the spectral (quasi-particle) density ^(k) is 
finite only at the four points in the Brillouin zone where the SC order parameter has nodes. 
In the striped superconductor it was shown that CDW order leads to an increase of spectral 
weight in the nodal direction while SDW order suppresses spectral weight in this region 
[65]. Here we observe that spectral weight reappears on most of the Fermi surface of the 
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Figure 2. Momentum-space quantities of the striped phase of a clean d-wave super- 
conductor, a) Momentum distribution n(k), (b) spectral density A(k), c) pair density P(k), and 
(d) spin density p,(k) for T = 0.025t,t' = -OAt.V = 1.5t,x = 1/8, U = 3.3t, and = 0. 

Horizontal and vertical axes are given in units of 7v/a. 

normal conducting state (Fig. 2 (b)). In the near nodal direction two peaks appear which 
are fingerprints of both CDW and SDW order with a dominance of the latter. Comparing 
to Ref. [65], one can attribute the peak around k (0.25,0.6)7r in the first quadrant of 
the Brillouin zone to spin order, while the peak Sit k (0.5,0.4)7r results from the charge 
order. These features bear resemblance to ARPES data found in La2-x-yNd^Sra;Cu04 and 
LSCO by Zhou et al [66]. Figure 2 (c) shows the pair density P(k) which accentuates the 
states contributing to superconductivity. The maxima in A(k) around k (0.5, 0.4)7r which 
originate from the CDW order [65] correspond to a suppression of the pairing in Fig. 2 (c). 
P(k) is maximal around k (±7r, 0) where the pairing amplitude is largest. These regions 
resemble those of a homogeneous dSC but stretching out less far into the nodal direction. 
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In contrast, P(k) is strongly suppressed at wavevectors k c::^ (0, ±7r) where the spin density 
is strong and exhibits local maxima around k c::^ (0,±3/47r). The unidirectional character 
of the spin-density wave is clearly seen in Ps(k) (see Fig. 2 (d)) which relates to the large 
amplitude oscillations of the staggered magnetization in real space, ^^(k) dominates regions 
in momentum space which P(k) would occupy in a non-magnetic homogeneous dSC. It 
extends far along the Fermi arcs, crossing the nodal points. Ps(k) has large intensities in 
the anti-nodal regions, i.e., it competes with superconductivity. Though maxima of Ps(k) 
are clearly separated from those of the pair density P(k), both quantities coexist in large 
parts of momentum space. 

2. Impurity-pinned stripes 

In the presence of a single strong impurity with potential strength V^^"^ = 10 1, stripes 
are pinned as is inferred from the order parameter pattern shown in Fig. 1 (b), where the 
impurity is located at the center. Impurity-pinned stripe solutions are robust against a 
change of the initial conditions, in contrast to the impurity-free stripe solution ((Fig. 1 (a)). 
Thus we expect fluctuating stripes in real materials to be pinned by impurities and to 
become static. Horizontally aligned stripes emerge as well in the magnetization and the 
electron density. The impurity pins a channel of reduced pairing- amplitude (Fig. 1 (b)), 
which minimizes the loss of pairing energy. On the other hand, inhomogeneities in the 
order parameter result in a charge density redistribution in the absence of particle-hole 
symmetry [58]. In fact, the channels of reduced pairing amplitude collect electrons, thereby 
shifting the electron density in these channels towards half-filling. This in turn favors the 
emergence of antiferromagnetic order in the regions of enhanced electron density. That 
is why the impurity pins one of the ridges of the staggered magnetization. Altogether, 
electron-rich stripes coincide with strongly magnetized stripes of reduced pairing amplitude. 
The characteristics of the impurity-pinned stripes are similar to those of the impurity-free 
case, except that the impurity-pinned stripes emerge independently of the initial conditions. 
Thus the momentum space quantities shown in Fig. 2 for the impurity-free stripe solution are 
essentially the same for the impurity-pinned stripes. Summarizing, we observe for sufficiently 
large on-site repulsion U > Uc-> as in the impurity- free case, the coexistence of a mdSC, a 
SDW, and a CDW with wavelength Aqdw — ^mdsc = 4a = 1/2Asdw pinned by a single 
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non-magnetic impurity. 

Real-space quantities of stripes pinned by a single strong non-magnetic impurity have 
already been investigated by Chen and Ting [67]. They obtained similar results for the 
[/-model at X = 15% hole doping which is assumed to be close to optimal doping. We also 
checked that the periodicity of the stripes for a single impurity at x = 1/10 remains the 
same as for x = 1/8. 

3. Vortex-pinned stripes 

The vortex-pinned stripes shown in Fig. 3 display essentially the same properties as 
those for the impurity- free and the impurity-pinned stripes. The amplitudes of the mdSC, 
CDW, and the SDW are marginally enhanced as compared to the impurity-pinned stripes 
because a vortex penetrating the superconductor acts as a far stronger perturbation than 
a single pointlike impurity. Just like the impurity, the vortices pin the ridges of the SDW. 
Stripes of reduced pairing run through the elliptically deformed vortices and thereby save 
condensation energy. The AF stripes are further attracted by the vortex since vortex cores in 
rf-wave superconductors strongly magnetize due to a spin-dependent splitting of the Andreev 
bound-state [38, 58]. Correspondingly, the vortex pins an electron-rich stripe with a filling 
shifted towards 1/2. From the Fourier transforms of the real-space quantities we obtain 
the same wavelength of the mdSC (Amdsc — 4a), the CDW (Aqdw — 4a), and the SDW 
(AsDW = 8a) at X = 1/10 and at x = 1/8. Similar results for coexisting CDWs and SDWs 
were observed in Ref. [63] within the [/-model for vortex-pinned stripes above a critical 
[/ at X = 15% hole doping. Hence the wavelengths of these unidirectional density waves 
are identical in the impurity-free and the impurity- and vortex-pinned stripes and robust 
against a change of hole doping or pairing interaction strength, as we verified by explicit 
calculations. 

4. Disorder-pinned stripes 

For the modeling of cuprates, such as BSCCO or those based on the parent compound 
La2Cu04, their dopant disorder has to be taken into account. In the following we consider 
explicitly intrinsic disorder and compare our results with experimental data in zero and finite 
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Figure 3. Vortex-pinned stripes. Results are shown for a magnetic flux $ = 2$o penetrating 
a 32a x 32a lattice which corresponds to a magnetic field B = 25.6 T. (a) SC order parameter, 
(b) electron density, (c) staggered magnetization (blue lines mark the zero-crossing). The model 
parameters were set to T = 0.025 t,f -OAt.V = 1.34t,x = 1/8, = 3.2 1. 

magnetic field. For modeling dopant disorder we set the impurity concentration equal to the 
hole doping = x. Moreover, we assume that the dopants' impurity potential decreases 
with increasing hole doping due to enhanced screening [53]. In Fig. 4 and Fig. 5 we show 
exemplary results for two specific impurity configurations which exhibit the characteristic 
behavior of the disordered stripe phase. We have verified that other impurity configurations 
yield density waves with the same wavelength as discussed here. For a Coloumb repulsion 
U = 3.3t we find quasi-unidirectional stripes at x = 1/8, slightly deformed by the presence 
of the non-magnetic impurities (see Fig. 4 (d)-(f)). Hole-rich paths (see Fig. 4 (e)) coincide 
with anti-phase domain walls, seen in the staggered magnetization profile in Fig. 4 (f), and 
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Figure 4. Emergence and stabilization of stripes in real-space quantities of d-wave 
superconductors by a shift of hole doping from 1/10 to 1/8 (nimp = x^U = 3.3 1). For 

this set of parameters density modulations emerge close to x = 1/10 doping (T/™p = 1.3t) in the 
presence of strong dopant disorder and become unidirectional by a change of hole doping towards 
1/8 = 0.9t). (a), (d) SC order parameter A^, (b), (e) electron density (c), (f) staggered 

magnetization (blue lines mark the zero-crossing) for T — 0.025 = —OAt^V — 1.6 1. 
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Figure 5. Fourier transformed SC order parameter (a), (c), electron density (b), (e), 
and magnetic structure factor (c), (f) for 1/10 (a)— (c) and 1/8 (d)— (f) hole doping. 

The dopants' impurity potential is T^mp = 0.9t at x = 1/8 and Vimp = 1.3 1 at x = 1/10. Note that 
the impurity sites are not considered in the FT electron density (b), (e) as described in Sec. II C. 
Parameters were fixed to ridop — x^U — 3.3t^T = 0.025 = —OAt^V = 1.6t. Horizontal and 
vertical axes are given in units of vr/a. 
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with stripes of strong superconductivity as observed in the SC order parameter Fig. 4 (d). 
Although the suppression of the electron density on the impurity sites dominates the charge 
pattern as seen in Fig. 4 (e), closer inspection reveals the existence of hole-rich channels 
(colored orange in Fig. 4 (e)), where (n^) is reduced compared to the average electron 
density. Additionally, the lines of zero staggered magnetization (blue lines in Fig. 4 (f)), 
which coincide with the maxima of the SC order parameter, also serve as a guide to the eye. 
The G?-wave order parameter (Fig. 4 (d)) varies in the strong-disorder limit with much larger 
amplitudes as compared to the previous cases, but it never changes sign. 

From the analysis of different disorder configurations we conclude that the hole-rich paths 
emerge only close to x = 1/8 hole doping and sharpen upon approaching x = 1/8. Simulta- 
neously the density modulations as a whole become quasi unidirectional. We find that, by 
shifting the hole doping towards 1/8, the system is driven into the stripe state. Moreover, 
we observe that the wavelength of the density modulation decreases by enhancing the hole 
doping from x = 1/10 to 1/8 (cf Fig. 5). This characteristic, which is also observed in 
neutron scattering experiments [68, 69], is not found in the impurity-free systems. Exper- 
imentally, SDWs with wavelengths Asdw(^ 1/10) 10a and Asdw(^ = 1/8) = 8a were 
inferred from the incommensurabilities. While this model calculation predicts exactly the 
same wavelength for the SDW at x = 1/8, the wavelength at x = 1/10 is slightly larger, 
that is AsDW — 12a, than in the experiment. Note that a decrease in hole density either re- 
duces the amplitudes of the density waves or enhances the wavelength. The latter is valid in 
this model calculation. Strongly disordered systems support the former, where the random 
impurity sites facilitate the setup of additional hole channels. Disorder allows the stripe 
pattern to adjust more easily to the average hole density, whereas in the clean case, the 
stripe pattern is tied to the underlying lattice. 

All real-space quantities verify horizontally (or vertically) oriented stripes at x = 1/8. 
The periodicity of these patterns is extracted from the Fourier transformed (FT) quantities 
as shown in Fig. 5. Intriguingly, the FT data in the strong disorder regime for x = 1/8 
(Fig. 5 (d)-(f)) peak in general at the same wavevectors as for the perfectly unidirectional 
stripes discussed for clean systems. The magnetic structure factor (Fig. 5 (f)) shows two 
dominating peaks at the incommensurate wavevectors = ^(1/2, 1/2 ± e) with incom- 
mensurability e = 1/8. The FT electron density (Fig. 5 (e)) and the FT SC order pa- 
rameter (Fig. 5 (d)) exhibit similar patterns with two broad maxima around wavevectors 
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Qc/p '::±2'K/a (0, ±5) with 5=1/4. The main difference to the clean stripe solutions is that 
the peaks in Fig. 5 (d), (e) are broadened, and a low-intensity substructure is observable. 
Obviously, the slight deviation from the characteristics of perfect unidirectional stripes is 
caused by the finite impurity concentration. We infer that coexisting mdSC, CDW, and 
SDW with wavelength Amdsc — -^cdw 4a and Asdw = 8a persist at x = 1/8 even in the 
regime of strong disorder. The wavelengths of the density modulations change with doping 
as is recognized by comparison of 1/10 and 1/8 hole doping. In Fig. 5 (a)-(c) the FT SC 
order parameter, the FT electron density, and the magnetic structure factor are displayed 
for X — 1/10. The maxima deviate clearly from those atx = l/8. At x = 1/10 the density 
modulations oscillate with wavelengths Asdw — 12a and Amdsc — 6a = Aqdw- In contrast 
to weakly disordered systems the incommensurabilities e and 5 in the strong disorder limit 
are sensitive to doping and increase with growing hole concentration. 

The momentum space quantities n(k), A(k), P(k), and Ps(k) (see Fig. 6) in the strong 
disorder limit resemble the results obtained for the impurity- free stripe solutions (see Fig. 2). 
Disorder slightly enhances the anisotropy, best visible in ^^(k). Importantly, in the disor- 
dered system continuous Fermi arcs appear in the near nodal direction (see Fig. 6 (b)). This 
agrees with ARPES measurements in LBCO [70], which has a strong dopant disorder. Due 
to the impurity potentials, through which the hole channels run, the CDW is stronger here 
as compared to the previously discussed cases, which results in the redistribution of spectral 
weight to the nodal directions [65]. 

Altogether, we find that a SDW, a CDW, and a mdSC coexist at x = 1/8 hole doping 
in strongly disordered systems such as dopant disordered LSCO. The density-modulations 
are pinned and stabilized, just as for a single impurity or a vortex, yet slightly deformed by 
the strong disorder. The doping dependence of the incommensurabilities e (x x and 5 oc x is 
similar to those observed in experiments [68, 69]. 

5. Density of states 

Characteristically all the identified stripe solutions exhibit a full gap in the density of 
states (DOS) (see Fig. 7). For comparison also the typical v-shaped gap of a homogeneous d- 
wave superconductor is shown in Fig. 7. In the presence of stripes a full gap opens around the 
Fermi energy. The full gap is a generic property of stripes in rf-wave superconductors which 
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Figure 6. Momentum-space quantities in the strong disorder limit, a) Momentum 
distribution n(k), (b) spectral density A(k), c) pair density P(k), and (d) spin density PsO<-) for 
T = 0.025 = -OAt.V = 1.6 t,U = 3.3t,x = 1/8 = ndop, 1^^'^''^ = 0.9 Horizontal and vertical 
axes are given in units of vr/a. 

persists also in the strong disorder limit. The full gap originates from a finite extended 
s-wave contribution in the striped rf-wave superconductor and the SDW gap present in 
the strongly magnetized stripes. In a homogeneous rf-wave superconductor the bond order 
parameter Aij changes sign between horizontal and vertical bonds which refiects the rf-wave 
symmetry of the energy gap Ak- Thus the summation of Aij over the vertical and horizontal 
bonds connecting to a given lattice site gives zero in a homogeneous rfSC. Due to the strong 
anisotropy of the stripe solutions, however, this summation yields a finite value here, which 
characterizes extended 5- wave superconductors. In a dSC coexisting with a SDW the gap 
is generically not centered around the Fermi energy, which is the reason why the full gap is 
asymmetric in all stripe solutions shown in Fig. 7. In addition, particle-hole symmetry is 
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Figure 7. DOS in d-wave superconductors in the presence of stripes. Results for the 
data sets used above for impurity-free, impurity-pinned, vortex-pinned, and disorder-stabilized 
stripe solutions for x = 1/8. The black curve corresponds to the DOS of a homogeneous d-wave 
superconductor. 

already broken by the constant random impurity potentials with l^mp > 0. Experimentally, 
particle-hole anisotropy was observed recently by angle-resolved photoemission spectroscopy 
(ARPES) in the pseudogap phase of Bi-2201 [6] which is attributed to competing orders 
and contrasted to homogeneous superconductivity. The appearance of a full gap was also 
found by Loder et al [41]. Obviously, the gap is largest in the unperturbed stripe solution 
{B = 0, Virnp = 0). This is due to the fact that with the impurity- and field-induced bound 
states, spectral weight is shifted into the energy gap. 

B. Impurities in the V-Model 

The l/-model presented by the Hamiltonian Eq. (6) was introduced in Ref. [41] as an 
alternative description for the coexistence of superconductivity with magnetic order. In the 
homogeneous case, the groundstate solution of the ]/-model for x = 1/8 is a mdSC state 
similar to Fig. 1(a), (c). The main difference is that the AF stripes are nearly maximally 
magnetized, i.e. the SDW oscillates with an amplitude ^ /xb, while Af goes to zero at the 
center of the AF stripes without changing sign. Thus magnetic order is dominating in the 
y-model, whereas superconductivity is reduced to straight, quasi one-dimensional lines on 
the domain walls between the AF stripes. This mdSC solution is degenerate with a "pure" 
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Figure 8. Stripe patterns in the disordered V-model. (a)-(c) Weak impurity potentials 
yimp ^ o.9t (x = 1/8 = ndop) and (d) strong impurity potentials = 2.0t (ndop = 4%). (a) d- 

wave order parameter A^. (b) Electron density, (c) Staggered magnetization, (d) Electron density 
exhibiting global phase seperation for strong impurity potentials. The turquoise dots indicate the 
positions of the impurities. Parameters were fixed to T = 0.01 1, = — 0.4 1, V = It. 

PDW, which features a periodic sign change in Af. Both, the "pure" PDW and the mdSC 
have identical pair densities P^(k) (c.f. Ref. [41]). 

In the strong disorder limit the U- and the l/-model show a severely different behavior. 
The solutions of the ]/-model can be divided into two regimes. If the scattering strength 
of the impurities is weak (i.e. 1/™p < t)^ the AF stripes remain straight, as shown in Fig. 8. 
The magnetic energy gain is dominant and prevents the superconducting stripes to wind 
around the impurities, which would allow to gain further condensation energy. Moreover, 
impurities in the straight one-dimensional SC channels strongly suppress superconductivity. 
This effect culminates in the disappearance of entire SC stripes, if the impurity density in 
their near vicinity becomes too large (Fig. 8 (a)). One of the neighboring AF domains (see 
Fig. 8 (c)) spreads over such a metallic line expelling the holes, which collect in the remaining 
superconducting stripes (Fig. 8 (b)). In this situation the AF domain- wall boundaries are 
shifted locally, according to the impurities' positions, but the AF spin stripes remain largely 
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straight (Fig. 8 (c)). The antiferromagnetic stripes are almost half filled and all sites within 
the AF stripes are nearly fully polarized, in contrast to the small amplitude modulations 
seen in Fig. 4 for the f/-model. Thus the superconducting stripes in the disordered V- 
model have a filling below 1/4 because the holes collect mainly in the remaining SC stripes. 
Similar to the SDW, also the CDW in the l/-model modulates with much larger amplitudes 
as compared to the t/-model. 

If the impurity potentials become larger than we expect that the impurity-pinning 
of stripes is energetically favorable as long as a mdSC state prevails, similar to the [/-model 
(c.f. Fig. 4(a)-(c)). However, we observe an overall disappearance of superconductivity 
accompanied by a global phase separation. All holes accumulate where the impurity con- 
centration is largest and antiferromagnetism weakest, while the surrounding region turns 
into a half-filled AF insulator (Fig. 8 (d)). The same occurs if weak impurities are numerous 
enough to destroy sufficiently many metallic (superconducting) lines so that the remaining 
ones cannot accommodate all the holes in the system. This behavior is characteristic for the 
strong-coupling limit described by the V-model: if the effectively one-dimensional striped 
system is forced to become two-dimensional because of disorder, the system divides into 
half-filled and empty regions featuring global phase separation. 

C. Discussion of the experimental observations in different cuprates 

We found that the [/-model allows for two-dimensional superconductivity in coexistence 
with a stripe order that adapts fiexibly to disorder, while the solutions of the l/-model remain 
quasi one-dimensional. Materials in which three-dimensional superconductivity coexists with 
stripe order, such as LBCO [7] or the rare-earth doped cuprate L8i2-x-y^dySixGn04 with 
X < 0.12 and y = 0.2 [1, 12], are certainly not as one-dimensional as the solution of the V- 
model. On the other hand the TZ-model reproduces the quasi one-dimensional characteristics 
of Nd-doped LSCO around y = 0.4, quite well [41], where the material becomes more 
anisotropic with increasing doping. As argued in Ref. [42], the models described by Hu 
T-Lv are the weak and strong coupling limits of a t-J like model, and the cuprates are most 
likely found somewhere in the intermediate coupling regime. 

The static, unidirectional CDWs and SDWs observed in LBCO and LNSCO [4] are often 
ascribed to the structural phase transition towards the LTT phase, where the LTT specific 

23 



buckling pattern of the CuOe octahedra induces an x-y anisotropy [71]. Our results imply 
that the strong dopant disorder, which is present in these substances, further stabilizes the 
stripe state. The results within the f/-model in the strong disorder limit reproduce typical 
properties of these substances well, such as the Fermi arc reconstruction in LBCO [70] 
or LSCO [72] and the doping dependence of the wavelength of the SDW as inferred from 
neutron scattering experiments in LSCO [68, 69]. 

Though spin stripes are observed in LSCO at x = 1/8, an ordering of the charges has so 
far no^ been detected by neutron scattering experiments [18, 19]. As LSCO does not exhibit 
an LTT phase we suggest that the experimentally observed spin stripes [18] are induced 
by the dopant disorder. The expected concomitant charge stripes are not found by neutron 
scattering experiments possibly because the charge modulations are weak and may therefore 
be concealed by the reduction of the electron density below the dopant sites. 

YBCO, which has neither an effective dopant disorder nor a LTT phase, exhibits dy- 
namic instead of static stripes [20, 73]. This agrees with the results of the [/-model for an 
impurity- free system as discussed. In the presence of a magnetic field charge stripes become 
static, while spin stripes are absent [23]. In this case the magnetic field induced charge 
stripes erase superconductivity. We find a similar behavior in the [/-model when we turn oflF 
superconductivity deliberately. The magnetization is greatly suppressed, while a dominant 
charge order survives. 

IV. SUMMARY 

The models we investigated in this article allow to explain a wide range of features of 
striped cuprate superconductors, which were found experimentally. Notably, the material 
specific characteristics of stripe phases in various cuprates can be modelled within the U- 
model. On the other hand, we found that AF correlations in the TZ-model tend to be stronger 
and the superconducting stripes are one-dimensional which reproduces some characteristics 
of Nd-doped LSCO quite well. 

Specifically, we found the coexistence of PDW's, SDW's, and CDW's with wavelengths 
AcDW — ApDW = 4a = 1/2 AsDW in ^ G?-wave superconductor for sufl&ciently large U in the U- 
model. Impurity-free, (single) impurity-pinned, and vortex-pinned straight stripe solutions 
exist over a broad doping range in the [/-model. Disorder dissolves the strict one-dimensional 
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stripe order into meandering stripe-like patterns. Only close to x = 1/8 hole doping, which 
is ideal for stripe formation, a unidirectional stripe order reemerges in strongly disordered 
systems. 

In the l/-model stripes are rigid and remain virtually undistorted in the presence of weak 
disorder. Superconducting stripes react very sensitively to disorder which may even re- 
move some of the superconducting channels. The remaining superconducting stripes exhibit 
strictly one-dimensional superconductivity appearing between straight AF stripes. Strong 
disorder, however, leads to global phase separation into dominantly antiferromagnetic re- 
gions with non-magnetic puddles absorbing the holes. 

The overall identical characteristics of the impurity-, vortex-pinned, and (strong) disorder- 
pinned stripes indicate that the stripes are not generated by impurities or fields. Instead, 
stripe states, which are energetically close to homogeneous solutions in unperturbed systems, 
are pinned and stabilized by inhomogeneities which are connected to charge redistributions. 

A generic feature of striped rf-wave superconductors in both the U- and the TZ-model 
is the opening of a full gap, induced by an extended 5-wave contribution. Without local 
probes for the superconducting order parameter it is not clear if pair-density modulations 
accompany the charge- and spin stripe order. . 
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